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The density of states of self-gravitational system diverges when the particles are spread to infinity. 
Other problem based an inhomogeneous distribution of particles,which motivate the gravitational 
interaction. In this sense the statistical mechanics of self-gravitational system is essentially an non- 
equilibrium problem. A new possible approach to statistical description of self-gravitational system 
has been proposed. The approach based on non-equilibrium statistical operator, which allow take 
into account inhomogeneous distribution particle and temperature in self gravitational system. The 
saddle point procedure, which used in the given method describes the spatially inhomogeneous 
distribution in self gravitational system accompanied by temperature changing. 
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The statistical description of a system of interacting particles is a difficult problem but has the permanent attention 
for the problems of astrophysics [jij.The self-gravitational system are interesting for testing ideas about the statistical 
mechanic description of systems governed by long range interaction. The statistical mechanics of self gravitational 
system turns out to be very different from that of other, more familiar, many-body systems. A self-gravitational system 
has also more general problem, studied for a long time ■ The standard methods of statistical mechanic cannot be 
carrier to study gravitational system. Due to this fundamental difference, the notion of equilibrium is not always 
well defined and those system exhibit a nontrivial behaviour with gravitational collapse. For system with long-range 
interaction the thermodynamically ensemble are inequivalent, negative specific heat Q in microcanonical ensemble 
which not exist in canonical description [4]. In self-gravitational system can increase entropy and equilibrium states are 
only local entropy maximum. However, if introduce a repulsive potential at short distance, complete core collapse is 
prevented and can to proved that a global entropy maximum now exist for all accessible values of energy Q ■ For the gas 
with pure gravitational interaction between its particles the partition function diverges, the energy of self-gravitational 
system are not as extensive parameter. In it was shown that a self-gravitating gas collapses. A nature of the collapse 
and its conditions are explained by using simple and clear consideration Q. The phase transition in such systems 
creates the problem of description in the mean-field thermodynamics approach Q. Two type approaches (statistical 
and thermodynamic) have been develop to determination the equilibrium states of self-gravitational system Q , [j| ■ 
About all this problem and possible solution very good described in review Q . The collapse in such system are begin 
as spatially homogeneous distribution of particles in the all system at once. Formation of the spatially inhomogeneous 
distribution of interaction particles is a typical problem in condensed matter physics and requires non-conventional 
methods of statistical description of the system was tailored to gravitational interacting particles with regard for 
an arbitrary spatially inhomogeneous particle distribution. Spatially inhomogeneous particle distribution can not 
describe in standard approach equilibrium statistical physics. The microcanonical ensemble is the correct description 
of an isolated Hamiltonian self gravitational system at fixed energy, the canonical ensemble is the correct description 
of a dissipative Brownian system at fixed temperature [6|]. For inhomogeneous distribution of particle can not fixed 
any energy or temperature. Temperature can change together with distribution of particle in particular if fixed 
energy of system. Very important question are in point can be considered inhomogeneous distribution of particle as 
equilibrium state. This state is not stationary and derivation of density produce flux of particle though changing 
of gravitational force which act on every particle. The statistical description of self-gravitational system must take 
into account the possible inhomogeneous distribution the temperature together with distribution of particles. A few 
model systems with interaction are known for which the partition function can be exactly evaluated, at least within 
thermodynamic limits 9] but not for inhomogeneous distribution of particle [l(|. Inhomogeneous distribution of 
particle and temperature need to use nonequilibrium approach to description of behaviour self gravitational system. 
In this article present the possible approach local equilibrium statistical operator [14J which leads to take into account 
the inhomogeneous distribution particle and temperature. Formation of the spatially inhomogeneous distribution of 



interaction particles requires a nonconventional method, such as use in , 12[ , [131 ] , which are based an additional 
field representation of statistical sum [151 ] . For progress the goal can use saddle point approximation which provided 
to nonlinear equation. Solution of this equation present real distribution particle and temperature under certain 
conditions. In this approach do not use the polytrophic dependence pressure from density. This dependence are result 
as consequence of the thermodynamic relation. Presented condition give the possibility determine statistical operator 
for self gravitational system and fully describe the thermodynamic behaviour of self-gravitational system. 
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Phenomenological thermodynamic based on the conservation lows for average value of physical parameter as number 
of particles, energy and impulse. Statistical thermodynamic nonequilibrium system based too on conservation lows not 
the average value dynamic variables but in particular for this dynamic variables. For determination thermodynamic 
function of nonequilibrium system are need use the presentation of corresponding statistical ensembles which take 
into account the nonequilibrium states of this systems. The conception of Gibbs ensembles can brings to description 
nonequilibrium stationary states of system. In this case can determine nonequilibrium ensemble as totality of system 
which be contained in same stationary external action. This system have same character of contact with thermostat 
and possess all possible value macroscopical parameters which compatibility present conditions. In system, which are in 
same stationary external condition will be formed local equilibrium stationary distribution. For exactly determination 



local equilibrium ensemble must accordingly determine the distribution function or statistical operator of system 14 1 . 

If assume that nonequilibrium states of system can determine through inhomogeneous distribution energy H(r) 
and number of particles (density n(r) = H r ~ r ) the local equilibrium distribution function for classical system 
can write in the form = Q z _1 exp{— J (/3(r)iJ(r) — rj(r)n(r)) dr} where statistical operator local equilibrium 

distribution can determine as: 

Ql = J DTexpS.- J \p{T)H{v)- n {T)n{v))dv\ (1) 

The integration in present formula must take over all phase space of system. Must note, that in the case local 
equilibrium distribution Lagrange multipliers /3(r) and 77(1") are function of spatial point. Can introduce the entropy 
local equilibrium distribution by ratio 

S = -Sp (/, In f l ) = \nQ, +J (P(r)H(r) - V (r)n(r)) dr (2) 

After determination of statistical operator can obtain all thermodynamic parameter nonequilibrium system. For this, 
can determine thermodynamic relation for the systems. The variation of statistical operator by Lagrange multipliers 
can write necessary thermodynamic relation in the form: 

- 6 -^ = {H(r)) l - V (r){n(r)) l (3) 

and 
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This relation is natural general prolongation, on the case inhomogeneous system, well-known relation which take 
place in the case equilibrium systems. The conservation number of particles and energy in system can present in form 
natural relations J n(r)dr = N and J H(r)dr = E. 

For further statistical description of nonequilibrium system is necessary determine Hamiltonian of system. In the 
case of self gravitational system Hamiltonian can present in the form: 

H(r) = ?^-n{v) + - f W(r,r')n(r)n(r')dr / (5) 



2m 2 _ 

where impulse density p(r) = "YlnPi5{Y — r{) and gravitation interaction energy can present in well-known form 

Cm 2 

G is gravitational constant and m is mass of particle. This Hamiltonian of system is possible to use if take into 
account consider moving in phase space not compressed gravitational fluid. It is valid for not collision systems and 
self gravitational system present obvious example of such system. 

In the case self gravitational system the nonequilibrium statistical operator can write in the form: 

Qi= J DT cxp^- J ^/3(r)^^ -r/(r)^ n(r)dr- ^ W(r, r')n(r)n(r')drdr'| (7) 
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In order to perform a formal integration in secondpart of this presentation, additional field variables can be introduced 
making use of the theory of Gaussian integrals 15 1, (l3j ]: 



expj-iy /3(r)W(r,r')n(r)n{r')drdr'\ = J D<pexp^~J /3(r)W _1 (r, r')<p(r)(p(r')drdr' - J y/j^<p(r)n(r)\ 

(8) 



where Dip 



UdVs 



and W 1 (r, r') is the inverse operator which satisfies the condition 



y/dct 2ir/3W(r,r') 

J dr'W - 1 (r, r')W(r', r") = 6(r — r"). The inverse operator M^ _1 (r,r') of the gravitational interaction, in contin- 
uum limit should be treated in the operator sense [it| . i.e. 



w _1 (ry) = — 



1 



r A r <5(r - r') 



(9) 



where A r - is Laplace operator in real space. After this manipulation the statistical operator can rewrite in the form 

-,2 



Qi 



DT / Dip exp 



" " (r) " n{r)dr ~ 8^G J (V ^ (r))2 dV } (10) 



In this functional integral now can be provide the integration on phase space. The integration over phase space can 
present as functional integration by DT 



Dn(r)Dp(r) 
(2ttK) 3 



with regard for the cell volume (27rfi,) 3 in the phase space of 

individual states 12| , and take into account the " Pauli principle" for classical system, that own spatial point can not oc- 
cupier two classical particle. As result can use the relation J exp(— A(r)n(r))Dn(r) — dct ^ r ) {l — exp (— J A(r)drj }. 
Now can make functional integration over impulse. If introduce new variable absolute chemical activity £(r) = exp rj(r) 
after integration over impulse can write statistical operator in the form: 



Qi 



Dip exp 



dr 



As shown before [ll|],[12| in all cases classical statistic £ < 1 and can use expansion ln (1 — £A) « —£,A- 
the nonequilibrium statistical operator in more simple form: 



Qi 



Dip exp 



^(V,(r)) 2 + e(r)^)lexp^W,(r) 



dr 



(11) 
and rewrite 

(12) 



In the case constant temperature /? and absolute chemical activity £ the statistical operator fully reconstruct the 
equilibrium grand canonical partition function [l2T | , 21 1 . 

The statistical operator allows obtain use the of efficient methods developed in the quantum field theory with- 
out imposing additional restrictions of integration over field variables or the perturbation theory. In our case the 
nonequilibrium statistical operator can rewrite in the form 



Qi = / D<pexp{-F(<p(r),S(r),0(r))} 



where effective nonequilibrium "free energy" can present as: 



F(¥>(r),£(r),/3(r)) = 



87rm 2 G 



(V^(r)) 2 +£(r) 



27rm 
h 2 (3(r) 



) ^ exp y/f3(7)ip(r) 



dr 



(13) 



(14) 



The functional F(ip(r), £(r), /3(r)) depends on distribution of the field variables ip(r, chemical activity £(r and tem- 
perature /3(r). The field variable contains the same information as original distribution, information about possible 
states of the systems. The saddle point method can now be further employed can find the asymptotic value of the 
statistical operator Qi; the dominant contribution is given by the states which satisfy the extreme condition for 
the functional. It's easy to see that saddle point equation present thermodynamic relation and it can write in the 
other form: as equation for field variable s ^ r ^ = 0, the normalization condition J rfr ^^ = — j g ^ ^ £(r))dr = N 



and low of conservation the energy of system — J dr s ^ T ^ ^(r) — E. Solution of this equation fully determine all 
thermodynamic parameter and present general solution for the behaviour of self gravitational system, whether this 
distribution of particles is spatially inhomogeneous or not. The spatially inhomogeneous solution of this equations 
corespondent the distribution of interacting particles. Very important note, that only in this approach can take into 
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account the inhomogeneous distribution of temperature, which can depend from spatial distribution of particle in 
system. In other approaches the dependence of temperature from spatial point was introduce trough polytrophic 
dependence temperature from density of particle in equation of state. In present approach this dependence leads from 
necessary thermodynamic condition and can be determine for different distribution of particles. From normalization 
condition J p(r)dr — N can introduce the density function by definition in the form p(r) = £ n^r) ) T ex P(v^W ( / 9 ( r )) 
and rewrite the necessary equation in more simple presentation. Equation for field variable can write in the form 



Atp(r) + r mV WK r ) = (15) 
and equation of conservation energy take the form: 

p(r) 



2(3(r) 



(3 - VW^(r))* = E (16) 



where r m = AirGra 2 . In the case absent any interaction, from normalization condition J ^^P±)^dV = N at constant 
temperature can obtain the the chemical activity £ = VWfi) - ^ an< l fr° m equation conservation energy relation 
between energy and temperature fiVT = E. If use obtained relation can determine usual entropy in the form 

S = In { ^ww ) 2 *lf ^a,t fully reproduce entropy micro canonical distribution gaseous phase for fixed number of 
particles and energy in system. In the case homogeneous distribution of particle of self gravitational system can take 
p = y — const. If temperature T is constant, the equation for field variable <p leads to Aip + 47r ^, m p = and have 
solution ip = ^y~-. Substitution this solution in the equation for conservation energy Jq (3 — ^-^)r 2 dr = E 

can obtain for the gravitational gaseous phase in limiting volume of size Rcan obtain that |iVT = E + ^g-. This 
relation completely reproduce well-known relation [6j which was obtained by standard approach. 

In general case the distribution of particle in self- gravitational system are inhomogeneous. Inhomogeneous dis- 
tribution of particle motivate the long-range gravitational interaction. Can consider other simple case. If suppose 
that exist polynomial relation between density and temperature as y/ /3(r) p(r) — A = const and density change in 
accordance with not divergence of number of particle, that possible if density descend as fourth time of distance from 
center p(r) — CY~ 4 from normalized condition can obtain, that C — NRq where Rq is smallest distance to center 
where classical particle create close packing structure. From equation by field variable (p(r) = that directly re- 
construct behavior of gravitational field. From this condition temperature descend as fourth time of distance and can 
write that yj j3(r) = or T = ^r~ s . From equation for energy can obtain relation between introducing constants 
in the form g^R? ~ ^fjp" = ^- In this simple case can obtain all necessary coefficient and spatial dependence density 
and temperature for inhomogeneous distribution particles of self gravitational system. 

Indeed, present nonequilibrium statistical description tell only possible dilute structure in self-gravitational system 
but not describe meta stable states and tell nothing about time scales a kinetic theory. The statistical operator have 
not any peculiarity for different value of gravitational field. The problem of description of the self-gravitational system 
of particles could be solved with current approach which take into account the inhomogeneous distribution of particles 
and temperature. For the first time turn out well described the formation spatial inhomogeneous distribution particle 
with accompanied by changing temperature such distribution of interacting particles. More over, the method used 
can be also applied for further development of physics of self- gravitational and similar systems. 



[1] W. C. Saslow, Gravitational physics of stellar and galactic system, Cambridge University Press, New York, (1987) 
[2] T. Padmanabhan, Statistical mechanics of gravitating system, Physical Report, 188, 285, (1990) 
[3] W. Thirring, Z. Phys. 235, 339, (1970). 

[4] Clement Sire, Pierre-Henri Chavanis, Phys.Rev.E. 66,046133,(2002) 
[5] Pierre-Henri Chavanis, Phys.Rev.E, 65,056123,(2002) 

[6] Pierre-Henri Chavanis, International Jourmal of Modern Physics B, 20, 3113, (2006) 
[7] E. B. Aronson and C. J. Hansen, Astrophysics Journal, 117, 145, (1972). 
[8] D. Lynden-Bell and R. Wood, Mon. Not. R. Astron. Soc. 138, 495 ,(1968). 
[9] R. Baxter, Exactly Solved Models in Statistical Mechanics, Acad. Press, New York, (1982) 
[10] J. J. Aly. Rhys. Rev. E.49,3771,(1994) 

[11] Y. D. Bilotsky and B. I. Lev, Teor. Math. Fiz, 60, 120,(1984) 

[12] B. I. Lev and A. Yu. Zhugaevich, Phys. Rev. E. 57, 6460 (1998). 

[13] H. Kleinert , Gauge Field in Condensed Matter, Word Scientific, Singapure,(1989). 



5 



[14] D. N. Zubarev, Non-equilibrium statistical thermodynamics (Consultans Bareu, New York), (1974) 

[15] R.L. Stratonovich, Sov. Phys. Dokl. 2, 416 (1984). 

[16] K. Huang , Statistical Mechanics, J.WileyandSons. New.York,(1963) 

[17] V. B. Magalinsky, JETP, 21, 167,(1965) 

[18] A. Isihara, Statistical Mechanics, State University of New York, (1971). 

[19] D. Ruelle, Statistical Mechanics. Rigorous Results, New York, Amsterdam, (1969). 

[20] M. Kiessling, J. Stat. Phys. 55, 203 (1989) 

[21] H. J. de Vega, N. Sanchez and F. Combes ,Phys. Rev.E, 54,6008,(1996), |arXiv:hep-th /9903236 v2, 26 Jun (2000) 
[22] A.D.Linde, Elementary particle physics and inflationary cosmology, Horwood Academic. Chan. Switzerland, (1990). 



